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,__( ' We have investigated the relation of the direction of the momentum among the matter, 

^ ' neutrino, and proto-neutron star in a collapse-driven supernova in order to discuss the pulsar 

, kick. In particular, we have investigated the effects of the pulsar motion on the explosion, 

■ which are neglected in the previous study. As a result, it is suggested that the direction of 
' the total momentum of the matter and neutrino is opposite to that of the momentum of the 

proto-neutron star in the asymmetric explosion models. This is because the center of the 
. explosion deviates from the center of the progenitor due to the pulsar motion. This picture 

^ ' is common among the asymmetric explosion models. So if we assume that the pulsar motion 

' is caused by an asymmetric supernova explosion, the neutron star born in SN 1987A, which 

, has not been found yet, will be moving in the southern part of the remnant. In other words, 

■ if we can find one neutron star in SN 1987A on the south part of the remnant, asymmetric 
T-H ' explosion models will be supported by the observation better than the binary models. 

^ ' 

Q> . §1. Introduction 

Q . It is a well known fact that pulsars in our galaxy have velocities much in excess 

\ of those of ordinary stars cP. It is reported that their transverse speeds range from 

^ ' to ~ 1500 km and their mean three-dimensional speed is 450 it 90 km s^^H*. 
• • . On the other hand, there are many theoretical models in order to explain the 

.5^ \ pulsar kick. One is that a neutron star in a binary system can escape _from the 

^ ' system with rapid speed due to a supernova explosion of the nascent staru). There 



. are also many mode 

into considerationQ-' 



) 



in which effects of asymmetric supernova explosion are taken 
For example, it is reported that a 1% anisotropic neutrino 
radiation from the proto-neutron star results in a kick velocity consistent with the 
observations!-*. However, there are few observations to determine which model is 
the most promising one. 

As for the asymmetric explosion models, only Burrows and Hayes (1996) per- 
formed numerical simulations and estimated the speed of the proto-neutron star in 
the supernova matter. They introduced the anisotropy of the system by artificially 
decreasing the density of the Chandrasekhar core within 20° of the pole, which may 
be realized by the convection during the stellar evolution. They also estimated the 
contributions due to the neutrino emission anisotropy and the ejecta motions. As 
a result, they reported that the proto-neutron star can be accelerated to ~ 530 km 
s~^. They also concluded that the direction of the total momentum of the matter 
and neutrino is opposite to that of the momentum of the proto-neutron star. 
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Recently, Nagataki (1999) discussed that the asymmetry of the observed hne 
profiles of Fe[II] in SN 1987A can be also explained by the asymmetry of the explo- 
sion. In the case of SN 1987A, more matter has to be conveyed to the north side 
than to the south side in order to reproduce the observed asymmetric line profiles. 
As a result, if we believe the result of Burrows and Hayes (1996), the proto-neutron 
star born in SN 1987A, which has not been found yet, will be moving in the southern 
part of the remnanti\ 

It is worth while discussing whether the conclusion of Burrows and Hayes (1996) 
on the momentum of the matter, neutrino, and pulsar is common among the asym- 
metric supernova explosion models. If so, the conclusion on the location of the 
neutron star born in SN 1987A will be valid as long as the asymmetric explosion 
models are believed. Moreover, if the neutron star is discovered on the south side as 
Nagataki (1999) predicted, the asymmetric explosion models will be supported by 
the observation than the binary models. 

In this paper, we investigate whether the conclusion derived by Burrows and 
Hayes (1996) on the momentum of the matter, neutrino, and pulsar is common 
among the asymmetric supernova explosion models. In particular, we investigate 
the effects of the pulsar motion on the explosion, which are neglected in the study of 
Burrows and Hayes (1996). In section 2, we show analytical estimates for the effects 
of the neutrino heating by the proto-neutron star. In section 3, numerical estimates 
of neutrino flux are shown. Discussion and conclusion are presented in section 4. 

§2. Analytical Estimates for the effects of the neutrino heating 

The mechanism of collapse-driven supernovae has been understood as follows 
: when the mass of the iron core of the progenitor exceeds the Chandrasekhar 
mass, the star begins to collapse. The collapse continues until the central density 
of the collapsing core reaches about (1.5-2) times the nuclear matter density {p = 
2.7 X lO^^g cm"'^), beyond which matter becomes too stiff to be compressed further. 
A shock wave then forms, propagates outward. At first, the shock wave is not so 
strong and stall at ~ 200 kmcl) in the iron core (it is called a stalled shock wave). 
However, by the continuous neutrino heating (~ 500 msH^), the shock wave is revived, 
begins to propagate outward again, and finally produces the supernova explosion. 
This phenomenon is called as the delayed explosion, which is the most promising 
theory for the mechanism of the collapse-driven supernova explosion. 

In this section, we investigate the contributions to the shock revival due to 
the injection of the momentum and the thermal energy by the neutrino heating, 
respectively. 

Behind the stalled shock wave, the equation of state (EOS) is well determined 
by the radiation pressure, the ideal gas pressure, and the pressure of degenerate 
electrons at zero temperature. To put it concretely, the pressure and the energy per 
unit volume are described as follows: 



p = aT'^/S + pkT/pniu + Pe{p) 



(2-1) 



Instruction for Making IMgX Compuscripts Using pT'pT]^ 



3 



pe = aT^ + pkT/{jg - l)^m„ + 3peip) (2-2) 

where a is the radiation constant, k Boltzmann constant, the atomic mass unit, 
fi the mean mass number of the gas, 7^ the adiabatic index of the gas. The pressure 
of degenerate electrons pe is given by 

Pe = 1.24 X lO^^dyn cm~^{Yep)'^^'^ (2-3) 

where Ye is the electron fraction of the systerrL J^he^ typical values for the p, T, p, 
jg, and Ye behind the stalled shock wave are!) 0)0) 10^° g cm'^, 1 MeV, 1, 1.2, 
0.4. The density dependence of the pressure is shown in Fig. 1. 
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Fig. 1. Density dependence of the pressure. Temperature and electron fraction are set to be 1 MeV 
and 0.45. Solid line: the radiation pressure. Dotted line: the ideal gas pressure. Short dashed 
line: the degenerate electron pressure. 



At first we consider the contribution of the injection of the thermal energy. From 
Eq. (2-1) and (2.2), the partial derivative of the pressure by the energy density can 
be written as follows: 



dp 
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(2-4) 



At the behind the stalled shock wave, the contribution of the radiation can be ne- 
glected (see Fig. 1). So we can rewrite Eq. (2-4) as 



dp 
Ye 



1 



{ig - ~ TP 



(2-5) 
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From Eq. (2-5), we can estimate the thermal pressure by the neutrino heating as 

p = pe/5. (2-6) 

Next, we consider the contribution of the injection of the momentum. That is, 
we have to estimate the ram pressure which is the added pressure when the flow is 
interrupted and the velocity becomes to be zero at the stalled shock front. However, 
we can not use the Bernoulli's equation behind the stalled shock wave. This is 
because the pressure is determined not only by the electron but also by the ideal gas 
at p = 10^'' g cm~^ (see Fig. 1). In other words, the pressure is not determined only 
by the density. So we have to give a rough order-estimation for the ram pressure, 
which will be enough for the discussion in this section. From the dimension analysis, 
we can estimate that the ram pressure is as follows: 

Pram ~ pv"^ (2-7) 

where Pram is the obtained ram pressure. 

The typical explosion energy of a supernova is ~ 10^^ ergsHH'. The total mass 
behind the stalled shock wave is estimated to be ~ O.IMq when the average density 
behind the stalled shock wave and the location of the stalled shock wave are assumed 
to be 10^*^ g cm~^ and 200 km, respectively. So the thermal energy density behind 
the stalled shock wave is estimated to be 5x10^^ erg g~^. So, using Eq. (2-6), the 
thermal pressure due to the neutrino heating is estimated to be 10^^ dyn cm~^. 
On the other hand, the velocity due to the momentum transfer from neutrinos to 
the matter is estimated to be 2x10^ cm s~^. Using Eq. (2-7), the ram pressure is 
estimated to be 4x10^^ dyn cm~^. 

In this section, we find that the contribution to the pressure due to the injection 
of the thermal energy will be greater than that due to the injection of the momentum 
by the neutrino heating. However, the ratio of the ram pressure relative to the 
thermal one (~ 1/25 in our analysis) is small relative to their absolute value. It will 
not be strange that their contributions become comparable in a realistic calculation. 
So we will investigate both of their effects on the explosion in the next section. 

§3. Numerical estimates of Neutrino flux 

3.1. Angular dependence of the normal component of the neutrino flux 

As we noted in introduction, the proto-neutron star might have moved from the 
center of the collapse due to the kick before the explosion. If the dynamics of the 
explosion is controlled by the thermal pressure of the matter heated by the emitted 
neutrino flux, the integrated flux determines the magnitude of the explosion. In this 
case it is expected that the neutron star (NS) heats the matter behind the stoled 
shock wave anisotropically and the irons would not be blown off spherically. In 
this subsection we show that there is an anisotropy of the integrated neutrino flux 
received on the back of the shock front and show how it depends on the azimuthal 
angel cos 9 to the line of the kick. 

For the simplification, we assume that the neutron star whose radius is = r„s 
moves at a constant velocity = Vns along the line of the kick from the center of the 
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collapse. In addition we assume that the neutron star starts to move as soon as it 
begins to emit the electron neutrino and it keeps radiating the neutrino at a constant 
flux To ( erg cm~^ sec~^). Here let t the passed time since the neutron star started 
to move. Then the distance from the center of NS to the shock front is given by 
r = Iq — Vnst, where is the radius of the shock front. The matter behind the stoled 
shock wave receives the neutrino flux ~ (rns/r)^ J-q per an unit time. Because NS 
has the rapid velocity enough to reach near the shock front until the explosion, we 
consider the integrated flux only for < t < tmax-, where tmax = (^o — Tns)/vns- 

Here we take the line of the kick as z-axis. As we show in Fig. 2, let r2 the 
distance from NS to an unit area behind the shock front and r] the angle between 
the direction to it and z-axis. Then the angle r] is expressed by 



V 




Fig. 2. Definition of coordinates and angles. The filled circle is the moving neutron star (NS). The 
outer circle represents the shock front. 



cos, = !:^Ml^^, (3-1) 

where r2 is given by ~ + (^o ~ ~ 2/o(^o ~ ^) cos6. If the neutron star keeps 
radiating the neutrino until it reaches the shock front, a net neutrino flux T±{cos 6) 
received on an unit area behind the shock front is expressed as the integrated normal 
component of the neutrino flux for < t < t^ax ■ Then it is given by 

2 



%ns WoV^L + 2/o(l-cos0)(/o-r„s)y ■ 
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If the neutron star does not move at all, the normal component of the integrated 
flux is equal anywhere behind the shock front and it is given by 



lo 



X tr. 



^ -^0(^0 - >'».s ')''ns 

Then the fraction of the received flux is estimated by 



v/rL + 2Zo(l-cos0)(Zo-r-„,)' 



(3-3) 



(3-4) 



In Fig. 3 we plot .?^x(cos 9)/ !F as& function of cos 9. Here we take the representative 
values r„5= 10km, Zo=200 km and t;„s=450 km sec~^. 
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Fig. 3. Plot of the integrated normal component of the neutrino flux received on the back of the 
shock front as a function of the azimuthal angle cos^. The integration of the flux is performed 
while NS is moving for < t < tmax- 

From Fig. 3 we find that if the explosion is triggered by the thermal energy of the 
matter heated by the neutrino, it is expected that for ^ > 7r/3 the matter is extremely 
blown off. Namely, a small amount of the iron is strongly emitted forward because 
of the large pressure gradient. On the other hand, the majority is mildly pushed 
backward because the gradient of the thermal pressure is smoother in the backward 
direction. In addition we can see that the asymmetry of the received flux is much 
more larger than 1% which is needed for the initial kick of NS. It ensures that we 
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can neglect the intrinsic anisotropy of the neutrino emission from the proto-neutron 
star as we assumed. 

3.2. Angular dependence of the forward component of the neutrino flux 

If the explosion is triggered by ram pressure caused by the momentum transfer 
from the emitted neutrino to the matter just behind the stoled shock wave, the 
magnitude of the explosion should be closely related to the forward component of 
the neutrino flux. In this subsection we show how the integrated forward (z-axis) 
component of the neutrino flux depends on the azimuthal angle cos 0. 

If the neutron star keeps radiating the neutrino until it reaches the shock front, 
the integrated forward component of the neutrino flux received on an unit area is 
given by 

^(cos 9)z = I J^o{ — ) cos rjdt 




r2 



ns 



+ 2^0(1 - COS 6*) 







(3-5) 



On the other hand, if the neutron star does not move at all, the integrated flux is 
given in Eq. (3-3) as before. Then the fraction of the integrated forward component 
is estimated by 

T,{cose) _ il ( 1 _i\ .3g. 

^ ^0 - Tns V ^/rls + 2/o(l - COS e)(lQ - rns) k)' 

In Fig. 4 we plot Tz{cos 6)/ as a function of cos 6. We find that if the explosion is 
caused by the momentum transfer from the neutrino flux to the matter behind the 
shock front, it blows off the matter forward only for 9 > tt/S and the majority of 
the matter is pushed backward, which is consistent with the conclusion derived in 
the previous subsection. This means that our conclusion on the flow of the matter 
is valid irrespective to the effects investigated. Also, because the asymmetry of the 
integrated neutrino flux is much larger than 1% almost everywhere, we should not 
be worried about the inherent asymmetric emission needed for the origin of kick. 
This means that the fact that the large amount of the matter is blown off backward 
does not depend on the mechanism of asymmetric supernova explosion. 

§4. Discussion and Conclusion 



We have discussed whether the conclusion derived by Burrows and Hayes (1996) 
on the momentum of the matter, neutrino, and pulsar is common among the asym- 
metric supernova explosion models. In particular, we have investigated the effects 
of the pulsar motion on the explosion, which arc neglected in their study. As a 
result, it was suggested by the present discussions that the direction of the total 
momentum of the matter and neutrino is opposite to that of the momentum of the 
proto-neutron star in the asymmetric explosion models. This is because the center 
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Fig. 4. Plot of the integrated forward component of the neutrino fiux received on the back of the 
shock front as a function of the azimuthal angle cos 9. The integration of the flux is performed 
while NS is moving for < t < t,nax- 



of the explosion deviates from the center of the progenitor due to the pulsar motion. 
This picture is common among the asymmetric explosion models. So if we assume 
that the pulsar motion is caused by an asymmetric supernova explosion, the neutron 
star in SN 1987A, which has not been found yet, will be moving in the southern part 
of the remnant. 

Moreover, the disccwery of the mystery spot, separated from SN 1987A to the 
south region by 60 masEj^, will support our conclusion. The existence of the bright 
spot on the south region will suggest that the existence of the strong shock wave on 
the south region, which is consistent with our conclusion (see Fig. 3 and 4). Here 
we have to comment on the second bright source which was reported to be detected 
on the north side of the remnantlilP. In their paper, it was concluded that the south 
spot has to be red-shifted and the north side spot has to be blue-shifted. It means 
that the south part of the ring around SN 1987A is nearer to us than the north part 
if these bright spots are ejected from the polar regions. This ccmcliision is opposite 
to the ones derived by the other many observations on the ringll^^ t3 . In this study, 
we have assumed that the north side of the ring is nearer to us and concluded that 
the neutron star born in SN 1987A is running in the south region of the remnant. 
However, if we believe the conclusion derived from the study of the second bright 
spot, our conclusion on the location of the neutron star is changed oppositely. We also 
give an additional comment on the mechanism of the jet-induced explosion which 
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may explain the existence of the mystery spotO tH^. These are the promising 
ones which may explain the existence of the mystery spot, but there has to be a 
mechanism which breaks the symmetry with respect to the equatorial plane. This is 
because the mystery spot seems to exist on only one side as discussed above (at least, 
the second bright source on the north side is much fainter than the one on the south 
side) . The effects of the pulsar motion on the asymmetric explosion discussed in this 
study explain such an asymmetry with respect to the equatorial plane naturally. 

As mentioned above, the asymmetric explosion models can explain naturally the 
asymmetry of the iron line profiles and the existence of the mystery spot at the same 
time. Also it is suggested that the pulsar born in SN 1987A is running in the south 
part of the remnant. On the other hand, the simple binary models in which the 
explosion is assumed to be spherically symmetric can not explain these observations 
and predicts that the pulsar is located at the center of the remnant. Moreover, the 
only possible candidate as a companion of the progenitor of SN 1987A is a compact 
object (neutron star or black hole) c3*. In this case, two compact objects will be found 
in the remnant of SN 1987A. One of them is found at the center of the remnant, if 
we believe the simple binary models. In other words, if we can find one neutron star 
in SN 1987A on the south part of the remnant and deny the existence of the another 
compact object at the center of the remnant, asymmetric explosion models will be 
supported by the observation better than the binary models. 

There may be a possibility to detect anisotropy in the young supernova ejecta 
using VLBA as long as it is located within 100 MpcElP. There may also be a 
possibility to find a pulsar in the ejecta. Increase of such observations will make it 
clear the relation between the asymmetric explosion and pulsar kick. 

In the delayed explosion model, about 1% of the neutrino energy is transfered 
into the energy of the explosion. In this situation we should delicately treat the 
neutrino transfer, and we should discuss such a subtle problem carefully. However, 
to qualitatively understand the physical mechanism in the model, we think that the 
neutrino flux and its integrated value give the adequate informations. Therefore in 
this study, we have investigated the effects of the pulsar motion on the explosion 
using simple analysis. It will be necessary to check our conclusion by performing 
realistic numerical simulations in which the effects of neutrino heating and its back 
reaction are taken into consideration. In the present circumstances, the effects of the 
pulsar motion in the iron core have not been taken into account and we have treated 
the neutrino flux as a simple thermal radiation in the numerical simulations con- 
cerning with the collapse-driven supernova explosion. However, we can understand 
its importance easily because the neutron star can reach to the stalled shock front 
in ~ 500 ms, which is the dynamical timescale of the stalled shock wave, as long 
as the neutron star moves with the observed mean speed (450 km s~^). Even if we 
assume that the pulsar is accelerated constantly and its velocity becomes to be 450 
km s~^ in 500 ms, the location of the pulsar at t = 500 ms becomes to be ~ 110 km 
from the center. In fact. Burrows and Hayes (1996) reported that the pulsar gets the 
velocity 500 km s"^ in 200 ms. We can easily guess that the effects of the pulsar 
motion on the dynamics of the explosion should be taken into consideration. Such 
calculations are now underway and we will report the results in the near future. 



10 K. Kohri and S. Nagataki 

Acknowledgements 

We would like to thank Dr. S. Yamada for useful discussions. This research 
has been supported in part by a Grant-in-Aid for the Center-of-Excellence (COE) 
Research (07CE2002) and for the Scientific Research Fund (199804502, 199908802) 
of the Ministry of Education, Science, Sports and Culture in Japan and by Japan 
Society for the Promotion of Science Postdoctoral Fellowships for Research Abroad. 

References 

[1] P.A. Harrison, A.G. Lyne, and B. Anderson, Men. Not. R. Astron. Soc. 261 (1993), 113. 

[2] A. Lyne and D.R. Lorimer, Nature 369 (1994), 127. 

[3] J.R. Gott, J.E. Gunn, and J.R Ostriker, Astrophys. J. 160 (1970), L91. 

[4] A. Burrows, and J. Hayes, Rhys. Rev. Lett. 76 (1996), 352 

[5] A. Kusenko, and G. Segre, Rhys. Rev. Lett. 77 (1996), 4872 

[6] S. Nagataki Astrophys. J. Suppl. 127 (2000), 1. 

[7] H.A. Bethe, Rev. of. Mod. Rhys. 62 (1990), 801 

[8] J.R. Wilson, in Numerical Astropysics, eds. D.G. Blair and M.J. Buckingham (World Sci- 
entific, Singapore). 

[9] A. Burrows, J. Hayes, and B.A. Fryxell, Astrophys. J. 450 (1995), 830. 



[10] C.L. Fryer, and A. Heger, |astro-ph/9907435 

[11] H.-T. Janka, E.M. Midler, Astrophys. J. 448 (1995), 109. 

[12] T. Shigeyama, K. Nomoto, and M. Hashimoto, Astron. & Astropys. 196 (1988), 141. 

[13] R. Nisenson, C. Rapaliolis, M. Karovska, and R. Noyes, Astrophys. J. 320 (1987), L15. 

[14] R. Nisenson, and C. Rapaliolis, Astrophys. J. 518 (1999), L19. 

[15] R.C. Rlait, R. Lundqvist, R.A. Chevalier, R.R. Kirshner, Astrophys. J. 439 (1995), 730. 

[16] E. Michael et al. Astrophys. J. 509 (1998), L117. 

[17] T. Riran, and T. Nakamura Nature 330 (1987), 28. 

[18] A.M. Khokholv et al. Astrophys. J. 524 (1999), 107. 

[19] A.L MacFadyen, and S.E. Woosley Astrophys. J. 524 (1999), 262. 

[20] T.J.B. Collins, A. Frank, J.E. Bjorkman, and M. Livio Astrophys. J. 512 (1999), 322. 



[21] B. Raczyhski, astro-ph/9909048 



